Curved Tails in Polymerization-Based Bacterial Motility 
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The curved actin "comet-tail" of the bacterium Listeria monocytogenes is a visually striking 
signature of actin polymerization-based motility. Similar actin tails are associated with Shigella 
flexneri, spotted-fever Rickettsiae, the Vaccinia virus, and vesicles and microspheres in related in 
vitro systems. We show that the torque required to produce the curvature in the tail can arise 
from randomly placed actin filaments pushing the bacterium or particle. We find that the curvature 
magnitude determines the number of actively pushing filaments, independent of viscosity and of the 
molecular details of force generation. The variation of the curvature with time can be used to infer 
the dynamics of actin filaments at the bacterial surface. 
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Introduction 

The bacteria L. monocytogenes, S. flexneri, the spot- 
ted fever group of Rickettsiae, and the Vaccinia virus are 
intracellular pathogens that move through the continual 
polymerization of actin ^, |[ ||] in distinctively curved 
"comet-tails" of actin filaments behind the motile parti- 
cles. While fascinating on its own, the actin comet-tail is 
functionally similar to the actin mesh in the lamellipodia 
of a locomoting eukaryotic cell, and the bacterial surface 
is analogous to the leading edge of the cell. Identifica- 
tion of the biochemical components involved has thus 
provided insight into the active regulation of actin poly- 
merization by the cell |^ — an essential cellular pro- 
cess 0. 

For polymerization-based motility, the force generated 
by actin polymerization at the moving object's surface 
drives the object forward against the viscous drag of the 
cytoplasm [g, |^. The necessary and sufficient bacte- 
rial contribution to motility is a single surface protein 
that orchestrates cellular cytoplasmic proteins to locally 
promote the nucleation, elongation, and cross-linking of 
actin filaments. In L. monocytogenes, this process is 
driven by the bacterial protein Act A [|o[ ^ 0, 
while 5*. flexneri expresses the protein IcsA for the same 
purpose Candidates for similar proteins have 

been proposed for spotted-fever Rickettsiae ^ and 
for the Vaccinia virus [T7| . Simplified systems that have 
been developed for the study of polymerization-based 
motility include Escherichia coli expressing IcsA on their 
surfacesjl^, |l^ and microspheres coated with purified 
Act A [|l9|. Similar motility mechanisms appear to be 
at work in endosomal rocketing [^0| , and in non-actin 
polymerization-based motility systems derived from ne- 
matode sperm pTf . Actin polymerization-based motility 
may even play an important role in vesicle trafficking 
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within the ceU ||2[ |3|. 

As the bacterium or particle is driven forward, a curved 
comet-like tail of actin filaments remains behind. Photo- 
bleaching experiments in L. monocytogenes and qual- 
itative observation of S. flexneri |l5| and of spotted-fever 
Rickettsiae |^ demonstrate that the tail is stationary 
with respect to the surrounding environment, probably 
due to steric or functional connections with the cellular 
cytoskeleton, so that the shape of the tail represents the 
path of the bacterium. The curvature of this path varies 
from bacterium to bacterium, and changes over time for 
individual bacteria. It is not known what determines the 
bacterial path, and hence the tail curvature, though no 
active control or chemotactic behavior has been proposed 
for these systems. 

Bacteria are functioning micromachines, but cannot be 
fully exploited without being fully understood. The L. 
monocytogenes motility system is well enough understood 
biochemically that ActA coated microspheres should 
reconstitute polymerization-based motility in solutions of 
purified proteins i.e. with total experimental con- 
trol. However, we do not yet quantitatively understand 
the motility enough to be able to use the bacterial or 
microsphere motion as a probe of the bacterial or cellu- 
lar conditions, or conversely to attempt to tailor those 
conditions to affect the bacterial motion. 

In this paper, we propose that the curvature results 
from the random location of actin filaments pushing 
against the bacterium. We show how the curvature of the 
bacterial path can be used to predict static and dynamic 
structure at the bacterial surface. The average curvature 
is determined by the number of active filaments pushing 
the bacterium. Information about filament lifetime and 
surface diffusion rates may be obtained from curvature 
autocorrelations, since curvature depends on the location 
of active filaments with respect to the bacterial surface. If 
filaments are closely localized to specific bacterial surface 
proteins, then lifetimes and diffusivities of those proteins 
can be inferred. We focus on the most common experi- 
mental geometry of a thin quasi-two-dimensional system 
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FIG. 1: The polar angle 9 and the axial angle with respect to 
the capsule-shaped bacterium. The bacterium is taken to be 
curving within the plane of the page. The half-length and the 
radius of the bacterium are a and 6, respectively. A schematic 
actin comet-tail is shown to the left, and the bacterium is 
pushed towards the right. We approximate the bacterium as 
a prolate ellipsoid for purposes of drag calculations. 



constrained between a glass slide and cover-slip, however 
we also discuss what would be expected for bulk geome- 
tries. In both cases, we discuss the apparent curvature 
appropriate to video microscopy. 



Curved Trajectories 

Curvature is defined as the rate of rotation of direction, 
K = d6/ds, or the rotation per unit path length, where 
6 is the polar angle in the current plane of motion (see 
Fig. 1^) and s is distance measured along the particle path. 
The radius of curvature, R = 1/\K\, equals the radius of 
the circle that locally best fits the path. Either K or R 
locally characterize the path, i.e. both may vary along 
the bacterial path. 

Curved trajectories imply a torque N is acting on the 
bacterium to balance viscous drag proportional to the 
angular rate of rotation, N — OCturm where Ctum = 
(327r/3)7y(a4 - 6'')/[(2a2 - b^)S - 2a] is a rotational drag 
coefficient for turning in 0, and 77 is the fluid viscosity, 
a and h are bacterial dimensions (see Fig. and S ~ 
2(a2 _ 62)-i/2 in[(a + ^[a^ - l?))lh\ [||. The viscous 
drag force F due to linear motion at speed v also carries 
a factor of the viscosity, where F = vfun and fu^ — 
16TrT]{a^-b^)/[{2a^-b^)S-2a] ||. Since the curvature 
of the path, K = 9/v = {N / Ctum) / (F / fun) , we can use 
the drag coefficients to obtain a remarkably simple result. 
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independent of the viscosity. This is fortunate since the 
effective viscosity of the cellular cytoplasm is strongly 
scale-dependent, ranging from 0.01 poise {10~^Pa s) for 
small loops on dye molecules p6| to 2100 poise (21QPa s) 
for 1.3^m diameter spheres [ |27| . We find that the curva- 
ture directly probes the ratio of force and torque applied 



to the bacterium, with an easily determined geometrical 
prefactor. 

The force, F — /qu, is proportional to the number 
of actively pushing filaments n, while the force per fila- 
ment /o depends on the specific details of the motility 
mechanism — as seen explicitly in thermal-ratchet mod- 
els of polymerization-based motility A complemen- 
tary coarse-grained elastic analysis of polymerization- 
based bacterial motion exists, however it is not conve- 
nient for determining curvatures. We take fo as constant 
in time, which amounts to considering only times much 
greater than the mean-time between actin monomer ad- 
dition or equivalently distances much greater than the 
monomer size 2.7nm. This is appropriate, since ob- 
served radii of curvature are larger than the bacterial 
scale (which is 1 — 2 microns). 

To calculate the torque N, we must consider the torque 
due to each filament. These individual torques depend on 
exactly where the filament pushes on the bacterium. We 
assume that the n actively pushing filaments are each 
randomly placed on the trailing end of the bacterium, 
so that each one will produce a random vectorial torque 
on the bacterium. The sum of many of these random 
torques will have a Gaussian distribution with zero mean. 
We can calculate the root-mean-square (RMS) torque, 
Nrms, from the local filament density. We take the grow- 
ing filament barbed-ends as uniformly distributed over 
the hemispherical cap at the end of the bacterium. In 
cross-section, this leads to enhanced filament density at 
the edges of the tail, similar to that seen in thin-section 
electron micrographs of L. monocytogenes ]2^ . This dis- 
tribution also follows naturally if the filament density fol- 
lows a uniform ActA surface density. In cross-sectional 
coordinates, where the cylindrical radius r ranges from 
to b, the filament density distribution is 
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where Pfrdrdcj) is the average number of filaments in the 
interval {r,r + dr) and {(j), (f> + dcj)). Here, a = n/(27r&^) is 
the uniform surface filament density on the hemispheri- 
cal end of the bacterium, and b is the bacterial radius. 
The mean-square torque perpendicular to the direction 
of motion is easily found 



d(f>rdr[{r sin < 



]Pfir, 



(3) 

where for sin and for cos are the two components of 
the torque, so that they add in quadrature. This leads to 
an RMS torque Nr^s = = fo9^JA^Ta/i. Taking 

the ratio Nrms/F to calculate the RMS curvature from 
Eq. |l|, we have 
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for a hemispherical distribution of filaments. Bacterial 
size and shape contribute to the curvature, as does the 
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average distribution of filaments Eq. ^ Other surface 
distributions of filament densities are also possible, and 
would affect the geometric prefactors in Eq. ^ though not 
the functional dependence on the number of filaments n. 

Remarkably, the average force per filament /o does 
not appear in our expression for the curvature, so that 
our results appear independent of the details of the force 
generation mechanism. However within polymerization 
ratchet models Q thermal fluctuations of semi-flexible 
actin filaments transverse to their length |^ could gen- 
erate transverse forces, which would lead to an increased 
torque and greater curvature than predicted a bove in 
Eqn. H This transverse contribution would depend on 
the biomechanics of the coupling between the actin fila- 
ment and the bacterial surface, which would also depend 
on the bacterial shape. Unfortunately it also depends on 
the effective and anisotropic elastic constants of the actin 
filaments | |29[ |, which in turn sensitively depends on how 
actively pushing actin filaments are cross-linked into the 
bacterial tail and the cytoskeleton — as can be see n by 
contrasting the elastic constants given by H and [|0|. 
Our simplified treatment corresponds to no coupling of 
forces transverse to the direction of bacterial motion. We 
hope that sensitive experiments can uncover the effects of 
these transverse forces and hence yield more insight into 
polymerization-based force generation, though we expect 
the effects to show up predominately in the geometrical 
prefactor or amplitude of the curvature in Eqn. ^ and 
not in the 1 / y/n dependence. We expect our subsequent 
analysis, on distributions and autocorrelations of the cur- 
vature, to be unaffected. 

The vectorial torque perpendicular to the bacterial di- 
rection of motion is Gaussian distributed, and the cur- 
vature is proportional to the magnitude of the torque. 
Within a bulk (3d) geometry the curvature has distribu- 
tion 

Pi{x) ^2xe-^\ (5) 
where x = \K\/Krms and Pj{x)dx = 1. 

Measuring curvature 

The curvature of the bacterial path, characterized by 
Eq. ^ and Eq. |^, is intrinsic to a given random placement 
of n filaments pushing against the bacterium. This intrin- 
sic curvature is constant for fixed filament locations on 
the bacterial surface. At any given time, the intrinsic cur- 
vature represents circular motion around the curvature 
axis. However, the instantaneous curvature axis is not 
necessarily parallel to the line of sight, so the apparent 
path would appear elliptical and have a non-uniform cur- 
vature. We must also consider dynamical effects which 
change the direction of the curvature axis, both through 
rotation of the bacterium and through motion of the fil- 
aments on the bacterial surface. 

Most experimental work to date has been done with 
restricted geometries, such as the typical gap of several 



microns |Q between a glass slide and its cover slip. For 
thin enough samples, curvature out of the plane will be 
restricted. If the axial angle measures the angle be- 
tween the vectorial torque and the normal to the sample 
plane in the line of sight, then the apparent curvature as 
measured from a microscope or video image will be 



Kapp = \K\ COS( 



(6) 



Positive and negative curvatures correspond to clockwise 
and counterclockwise curved paths in the microscope im- 
age, respectively. 

A single bacterium with a fixed intrinsic curvature, 
would eventually uniformly explore </> £ [0, 27r] through 
rotational diffusion As it does, the apparent curva- 
ture will change. For a fixed intrinsic curvature, sampling 
at uniform time intervals, we would measure a distribu- 
tion of apparent curvatures given by 



(7) 



where a:= \Kapp\/\K\ e [0,1]. 



If an individual bacterium changes its intrinsic curva- 
ture in time, then over sufficient time K will explore the 
entire intrinsic distribution, Eq. |[ The ensemble distri- 
bution Pe{x) of apparent curvatures will then be given 
by Eq. convoluted with Eq. ||: 



Pe{x) 
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where x = \Kapp\/ {K.rms/ "^)^ where the rms-apparent 
torque is equal to Krmsl^f^- Pe{x) also follows directly 
from the Gaussian distribution of each component of the 
torque. This "ensemble" distribution also characterizes 
the apparent curvatures of large groups of bacteria, since 
they will each have a different intrinsic curvature chosen 
from Eq. |^ and each will have a random axial angle (/). 
The differences between the ensemble, apparent, and in- 
trinsic curvature distributions is dramatic, as shown in 
Fig.|. 

For bacteria in a bulk sample, two angles are needed to 
characterize the bacterial path with respect to the viewer 
— a polar angle 5* with respect to the line of sight, and 
an azimuthal angle about the line of sight. One must 
also specify the axial angle about the bacterial axis. If 
we line up the bacterium to travel away from us along 
the line of sight, before rotating it towards us by ^ in 
the (po azimuthal direction, and measure </) as the angle 
between the vectorial torque and (j)o, then the apparent 
curvature is 



app,3d 



~ ~\K\ cos( 



'sin^ ^. 



(9) 



This is the curvature as measured from a microscope or 
video image, and takes values from (—00,00) with no 
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FIG. 2: The curvature distributions Pe{x) [thick solid], Pa{x) 
[solid], and Pi{x) [dashed] for ensemble, apparent, and intrin- 
sic curvatures, respectively, vs the scaled curvature, x (see 
text). We expect Pe(x) if a single bacterium is tracked over 
large times or if a collection of bacteria are tracked, Pa{x) 
if enough data can be gathered for a single bacterium before 
the intrinsic curvature changes, and Pi{x) for a single bac- 
terium if the curvature is extracted from a sufficiently long 
3d trajectory in a bulk sample. 



explicit 4>o dependence. Extremely large apparent curva- 
tures are seen when bacteria are moving towards or away 
from the observer with correspondingly small velocities, 
with small. The apparent speed is Vapp = f sin^f, 
where v is the bacterial speed along its path. A conve- 
nient quantity is obtained by "normalizing" the apparent 
curvature by multiplying by v'^ , to obtain 



app — Kapp^M'^app 



(10) 



where K = \Kv'^\. This also decreases the weight placed 
on apparently stationary bacteria, which can be hard to 
distinguish from other objects, and simplifies the analy- 
sis. Kapp will have the ensemble distribution Pe{x) if (p 
is uniformly explored in time. Note that time or path- 
length weighted distributions will differ in bulk samples, 
since the apparent speed Vapp will vary dramatically even 
for a constant intrinsic bacterial speed. The distributions 
presented in this paper are for time-weighted sampling, 
appropriate for video microscopy and/or for planar (2c?) 
geometries. 



Curvature Dynamics 

For individual bacteria, the experimentally apparent 
variation of curvature from one moment to the next is 
striking [|[ Restricting ourselves to thin planar 

samples, the angular diffusion of the bacterial orienta- 
tion can lead to changing apparent curvatures through 
changing axial angles 4> in Eq. |l^, and the intrinsic cur- 
vature can also vary if filament locations on the bacte- 
rial surface move significantly over time. We explore two 
cases, where active (pushing) filaments are removed and 



randomly replaced on the bacterial surface, and where 
active filament locations randomly diffuse. We will later 
apply these cases to characterize filament repositioning 
for different motile systems. 

It is easiest to characterize changes in the net torque 
acting on the bacterium with respect to a reference frame 
fixed to the bacterium. We consider the correlation of the 
intrinsic vector torques separated in time by At > 0: 



AN{At) = (Nit) -Nit + At)) 



-At/T 



(11) 



where a static intrinsic curvature corresponds to t = oo, 
and the average is over the initial time t. The second 
line follows directly if each filament has a lifetime r after 
which it is replaced randomly on the rear of the bac- 
terium by another filament. If new filaments are ran- 
domly placed they will be uncorrelated with other fila- 
ments. The correlation will then be proportional to the 
fraction of filaments that have not been replaced between 
the two times, i.e. e"^*^"*^^/'^. Exponential decay also 
applies for actin filaments whose fast-growing barbed-end 
positions diffuse over the bacterial surface with diffusion 
constant D. Solutions of diffusion on a spherical sur- 
face of radius b leads to t = b^/{2D). These results 
would apply directly to proteolysis/replacement or sur- 
face motion of bacterial proteins such as ActA or IcsA 
if active filament positions are localized to such bacterial 
surface features (see below). They should also apply to 
dynamics intrinsic to the actin tail through capping and 
nucleation of active actin filaments , where capping is 
a loss mechanism. In this case, the timescale of of auto- 
correlation decay would depend on the details of filament 
nucleation — how new filaments are placed with respect 
to pre-existing active filaments. 

Azimuthal diffusion will not affect the intrinsic curva- 
ture or its correlations, but it can contribute to decay 
of correlations of the apparent curvature by changing the 
apparent curvature over time. Azimuthal diffusion obeys 
((A(/))^) — 2DA^t, where Acf) is the net angle of rotation 
in the time At and Da is the diffusion constant. This 
has a direct effect on the correlation between apparent 
curvatures separated in time by At > 0: 

AappiAt) = {Kapp{t)Kapp{t + At)) 

= (i^2^,)e-^*/^(cos0(i) cosm) + A0)) 
= (i^L.)e-^*/^e-^-^V2, (12) 



where r is the decay time from intrinsic correlations in 
Eq. ^ and the average is over t for a single bacterium. 
[We have used the identity (e") = e^^^ for Gaus- 
sian distributed x.] Of course, rotation of the bacterium 
around its long axis will lead to decaying autocorrela- 
tions only if active filament tips are localized to bacterial 
surface features. If not, azimuthal diffusion might not 
affect the apparent curvature. 

For individual bacteria tracked for times much less 
than r, the intrinsic curvature will appear constant. For 
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times much longer than r, each bacterium will sample 
the ensemble of intrinsic curvatures. The characteristic 
timescale may be measured from the decay of curvature 
autocorrelations. Filament decay/replacement, filament 
diffusion, and axial rotation all contribute to exponential 
decay of the apparent curvature correlations. Their con- 
tributions to a particular motile system may be separated 
through independent measurements or through system- 
atic studies where parameters such as the particle size or 
the cytoplasmic viscosity are varied. 

Random rotation and diffusion of bacterial positions 
will also contribute to the measured curvature and cur- 
vature autocorrelations. In principle this is a compli- 
cated hydrodynamic effect leading to exponential 
asymptotic decay, however the timescales are very short 
(r ^ R^p/f] ^ 10~^s for a bacterium in a cell, where R 
is the cell size and p the cytoplasmic density) compared 
to the measurement intervals in typical bacterial exper- 
iments (seconds). The autocorrelation decay will effec- 
tively be discontinuous at At — 0, where thermal and 
measurement jitter will contribute at At — but not for 
At > 0. To eliminate those contributions, the experi- 
mental RMS curvature Krms should be extracted from 
the At limit of the autocorrelations or should be 

fit from sufficiently long segments of the bacterial path. 

Curvature autocorrelations in bulk samples are simple 
only when the intrinsic curvature is extracted from full 
3d tracking of the bacterial trajectory (see e.g. @). In 
that case, Eq. O will describe the autocorrelation decay. 



Motile systems 

In this section we discuss several specific motile sys- 
tems, and use the details to refine our discussion of curved 
bacterial paths. For illustrative purposes, curvature has 
been estimated from published images of L. monocyto- 
genes. This should be considered an order of magnitude 
estimate only. We only analyze motion within Xeno- 
pus laevis cell extracts, since cellular organelles and cell 
membranes, which can locally affect bacterial trajecto- 
ries through collisions and which are hard to control for 
in published images, are absent. Also absent in extracts 
is a polarized cytoskeleton, which could plausibly align 
bacterial motion in intact cells — this could be explored 
through a systematic comparison of bacterial motion in 
cells and in cell extracts. It must be emphasized that 
proper studies of curvature require unbiased data and in- 
dividual images previously selected for publication may 
be biased by aesthetic considerations. Distributions and 
autocorrelations require much more data than is avail- 
able from published individual images, and will require 
analysis of video data. 

The details of the nanoscale mechanical connection be- 
tween the actin filament tail and a particular bacterium 
or motile particle are not yet known, nor are the details 
of the dynamics. Indeed, these details may differ for dif- 
ferent bacteria or for different natural or reconstituted 



cytoplasmic environments. We present some plausible 
scenarios below and indicate the expected results of a 
curvature analysis in each. 



Listeria monocytogenes 

Motile L. monocytogenes have a distribution of tail 
lengths ranging up to about 15 pm ||2^, and speeds 
of up to 0.4 pm/s p^ . Mature L. monocytogenes are 
roughly cylindrical Gram-positive bacteria, 1.5 pm long 
with a diameter of approximately 0.5pm |32|] . Polar sur- 
face expression of ActA is required llfl, l lll for motility. 
L. monocytogenes in Xenopus extract ||37|| liave apparent 
curvatures of approximately K « 0.08^m~^, correspond- 
ing to rt = 20 filaments pushing on the bacterium. While 
this is a relatively small number, it is consistent with 
electron microscopy images . 

The role of ActA in polymerization is being uncovered 
psf but questions remain. It is not yet clear whether 
individual filaments are associated with individual ActA 
molecules and if so, for how long. Direct mechanical 
attachment is possible and indeed indicated by optical 
tweezer studies Indirect attachment is also possi- 
ble if the complex of ActA and cytoplasmic proteins 
serve as a source of monomeric actin with locally en- 
hanced polymerization affinity. For example, profilin has 
been shown to interact with ActA through vasodilator- 
stimulated phosphoprotein (VASP) ||lf| and hence can 
provide a local "plume" of profilin- ATP-G-actin which in 
some conditions can polymerize more readily than ATP- 
G-actin . We can readily obtain the leading behavior 
of the steady-state concentration C diffusing a distance 
d from a disk-like source of radius s and strength Co [0 , 
C/Co = 2s /{nd). The diffusive plume would provide sig- 
nificantly enhanced polymerization only to distances on 
order the size s of the ActA itself. If ActA is well sep- 
arated on the bacterial surface, then filament nucleation 
and growth would be closely associated with individual 
ActA molecules, which in turn would be held station- 
ary by the external peptidoglycan layer of the bacterium. 
The proteolysis induced lifetime r of individual ActA is 
approximately 2 hours in vivo [Q, and would directly 
contribute to the curvature autocorrelation decay given 
by Eq. 

If there is a connection between the bacterium and its 
tail j^], then axial diffusion will be dramatically reduced 
and autocorrelation decay due to proteolysis should dom- 
inate and could be directly observed. However, if the 
bacterium is not tightly attached to its tail, then it will 
rotate with Da — kBT/CA ~ lQ~^rad^ s~^ , correspond- 
ing to a timescale in Eq. Oof 1/ Da ~ 1000s ^ r, where 
Ca = (327r/3)77(a2 - b^)^{2a - b'^S) is the viscous drag 
coefficient for axial rotation and we take 77 = 30 
poise (3Pa s) following [||. 
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Shigella flexneri 

The Gram-negative S. flexneri has a similar motiUty 
mechanism to L. monocytogenes [p^ |4^ ; for example, 
it has a unipolar surface protein required for motility, 
IcsA. S. flexneri are about 2.3^m long and 0.5/iTO in di- 
ameter, and move at speeds comparable to L. monocyto- 
genes |3^. However, differences are observed between S. 
flexneri and L. monocytogenes. The tails of S. flexneri 
appear to have fewer actin filaments than L. monocyto- 
genes |4^. IcsA is targeted to one bacterial pole in 5*. 
flexneri and may diffuse in the outer membrane |44 , in 



contrast to L. monocytogenes where ActA is stationary. 

Curvature studies can help investigate these differ- 
ences. For example, if fewer filaments are actually push- 
ing the bacterium, rather than simply fewer filaments 
involved in cross-linking, then according to Eq. ^ the 
curvature of the path of 5". flexneri will be systemati- 
cally larger. If actin filament tips are associated with 
individual IcsA, then diffusion of IcsA on the outer bac- 
terial membrane will contribute a term e"^*/'^'*'^-' where 
Tdiff = b'^/{2D) in Eq. |l^, in addition to the finite IcsA 
lifetime due to proteolysis [Q . Unfortunately, S. flexneri 
are not motile in Xenopus extracts jl^, ^ , but they do 
have qualitatively similar curvature to L. monocytogenes 
in intact cells fB2|. 



Spotted-fever Rickettsiae 

The spotted-fever group of Rickettsiae use actin-based 
motility for intracellular movement ||l6|] . While a surface 
protein (rOmpA) of motile R. rickettsii has been impli- 
cated in tail formation and has sequence similarity to a 
domain of IcsA |Q , its surface distribution and specific 
biochemical role have not yet been characterized. The 
most studied species, R. conorii Q and R. rickettsii 
are roughly the same size as L. monocytogenes, but move 
only 1/3 as fast. 

While Rickettsiae are obligate pathogens that are not 
viable in Xenopus extracts, the curvature of their paths in 
intact cells is qualitatively smaller than that of L. mono- 
cytogenes or S. flexneri |Q. Yet the tails in R. conorii 
are found to have very long parallel filaments, with rel- 
atively few filaments observed in cross-section electron- 
micrographs ]3^ . It is worth considering two possibil- 
ities. The first is that the filaments are not uniformly 
distributed on the bacterial surface. The numeric pref- 
actor of the curvature Eq. ^ will range from for polar 
filaments, to •\/3/2 for uniform distribution on a hemi- 
sphere (Eq. H), to 3/2 if all filaments are along the outer 
edge at r = 6. A smaller than expected curvature can 
result from clustering of active filaments near the pole. 
However this appears to be unlikely from the electron 

A second pos- 



micrographs of decorated actin tails |32 



sibility is that the filaments are not randomly distributed 
on the bacterial surface. If the arrangement is symmet- 
ric, or regularly spaced, then the expected torque would 



be reduced. A paracrystalline surface structure is ob- 
served on Rickettsiae bacteria , though it is unknown 
whether it affects polymerization dynamics. 



Viruses and Vesicles 

Actin polymerization-based motility is exhibited by 
the Vaccinia virus, where ~ 200nm diameter oblate virus 
particles move at 0.05fim/s with tails of length 8/im in 
HeLa cells |^ . One of the interesting unsolved puzzles of 
this system is that the virus always travels in the sym- 
metry direction, which is the orientation of highest drag. 

If filaments are localized to the viral surface by a 
specific protein and if it diffuses over the viral 

surface, the intrinsic curvature will change with time 
with T = b'^/{2D) following Eq. |ll| and iQ, where b 
is the vesicular radius and D is the diffusion constant 
of the motility protein on the viral membrane. Using 
D sa lO^^cm^/s, appropriate for diffusion within vesicu- 
lar membranes, we have t = 0.05s! This is much less than 
the timescale of rotational diffusion and would dominate 
autocorrelation decay if present. 

Polymerization-based motility with curved comet-tails 
also occurs in motile vesicle systems, which are attractive 
systems for systematic study since they have fluid outer 
layers, spherical geometry, and a variety of sizes. Vesicle 
motion has been re-constituted in Xenopus egg extracts 
p2| , in endocytosed vesicles , and in extracts of ne- 
matode sperm |2^. In some systems the vesicle lipids 
directly mediate actin polymerization [^2[ |2^, in which 
case the filaments are unlikely to be localized to particu- 
lar spots on the vesicle surface. Effective filament motion 
could still occur due to random nucleation and loss of fil- 
aments from the vesicle surface, and this is a possible 
mechanism of intrinsic curvature autocorrelation decay 
in bacterial systems as well. 



Microspheres 

A simplified in vitro system using small polystyrene 
microspheres, coated with purified ActA and added to 
Xenopus lavis egg extract, has been shown to reconsti- 
tute actin-based motility jl^. Could different intrinsic 
curvatures be sampled by a single "inert" microsphere? 
The ActA has its transmembrane domain replaced by a 
6x His repeat, and is non-specifically bound to the car- 
boxylated polystyrene microsphere. It may be possible 
for the ActA to randomly crawl on the microsphere sur- 
face without detaching (for example, see Allowing 
a small surface diffusion rate, D = lO^^^cm^/s, results 
in a decay time for correlations of r = b'^/{2D) ~ 325s 
for a 0.5/zm diameter microsphere. Data from curvature 
studies may thus may be able to demonstrate diffusion 
of ActA on the microsphere. 

Microspheres are also good systems for systematic 
studies of various radii b with a constant surface density 
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of ActA and, presumably, actin filaments cr. In the au- 
tocorrelation decay, a finite filament lifetime would make 
a contribution that scales as t ^ const., diffusion on 
the microsphere surface would have r ~ 6^, while axial 
diffusion would have 1/Da ~ |2^- Varying the cyto- 
plasmic viscosity, on the other hand, should only affect 
axial diffusion, with \/Da ^ 1]- 

Summary 

Random filament interaction with the bacterial or par- 
ticle surface can explain the characteristic curved paths 
observed in polymerization-based motility systems, such 
as L. monocytogenes, S. flexneri, spotted-fever Rick- 
ettsiae, Vaccinia virus, and motile lipid vesicles and mi- 
crospheres. We distinguish between the intrinsic curva- 
ture, which can only be measured with the full three- 
dimensional trajectory of the bacterium or particle, and 
the apparent curvature observed in microscope images. 
We derived explicit distributions for these curvatures and 
showed how they can uncover important qualitative dif- 
ferences between the various polymerization-based motil- 
ity systems. We showed, in Eq. |2|, how the lifetime and 
dynamics of surface-associated proteins, such as ActA or 
IcsA, affects the evolution of the intrinsic curvature of the 
motion of individual bacteria, virus particles, vesicles, or 
protein-coated microspheres. 

Systematic experimental analysis of curvature in 
polymerization-based motility systems has not yet been 
done, but would supplement genetics, biochemistry, and 
microscopy by providing structural information about 
the interface between the actin tail and the bacterium. 
Curvature studies could estimate the number of filaments 
actively pushing the bacterium, the distribution of these 
active filaments on the bacterial surface, their localiza- 
tion with respect to motility protein complexes, and sur- 
face motility protein lifetime and diffusion on the bac- 
terial surface. A similar analysis can be done for virus, 
vesicle, and microsphere systems. 

It is useful to summarize the specific applications of 
this analysis. There are three. First, the relation be- 
tween curvature and number of filaments in the actin 
tail, Eqn. ^, can be compared with electron-microscopy 
cross-sections and with normalized fluorescence studies. 
The relation between curvature and particle size, also in 



Eqn. 1^, can be used for different sized particles with sim- 
ilar surface preparations, such as microspheres. Second, 
the distribution of observed curvatures is predicted to be 
Gaussian, Eqn. ^. For particles with constant intrinsic 
curvatures a qualitatively distinct distribution of appar- 
ent curvatures, Eqn. 0, is expected. If three-dimensional 
tracking of particles in a thick sample is used, then the in- 
trinsic distribution would apply, Eqn. |^. Third, the vari- 
ation of curvature in time is predicted to be described 
by an exponentially decaying autocorrelation function, 
Eqn. |T^. The timescale of autocorrelation decay, t, char- 
acterizes how the position of active filaments in the bac- 
terial tail change. Studies of different size particles, or 
direct tracking of azimuthal particle rotation, could help 
to untangle the possible contributions to curvature auto- 
correlation decay. 

We have made simplifying assumptions to facilitate our 
analysis. We have assumed particular particle shapes and 
surface distributions of filaments. We also took the indi- 
vidual filament forces /o to be in the direction of particle 
travel. Different shapes, surface distributions, and fila- 
ment orientations would change the numerical prefactor 
in Eqn. ^ though the curvature distributions would not 
be affected. We have also assumed that individual fila- 
ments are independently randomly located. Non-random 
symmetric filament locations will result in curvatures 
much less than predicted in this paper. In contrast fila- 
ment distributions that not symmetric will lead to qual- 
itatively larger curvature than discussed here, where we 
assume a random but azimuthally symmetric distribu- 
tion. We have also assumed that the viscosity does not 
vary strongly over bacterial length-scales. Strong viscous 
heterogeneities, perhaps caused by particle motion itself, 
as well as local constraints posed by cellular organelles 
and membranes, will affect particle trajectories. This 
could dominate the effects described here in vivo. 
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